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Abstract — Fading MIMO relay channels are studied analyti- 
cally, when the source and destination are equipped with multiple 
antennas and the relays have a single one. Compact closed-form 
expressions are obtained for the outage probability under i.i.d. 
and correlated Rayleigh-fading links. Low-outage approxima- 
tions are derived, which reveal a number of insights, including the 
impact of correlation, of the number of antennas, of relay noise 
and of relaying protocol. The effect of correlation is shown to 
be negligible, unless the channel becomes almost fully correlated. 
The SNR loss of relay fading channels compared to the AWGN 
channel is quantified. The SNR-asymptotic diversity-multiplexing 
tradeoff (DMT) is obtained for a broad class of fading distri- 
butions, including, as special cases, Rayleigh, Rice, Nakagami, 
WeibuU, which may be non-identical, spatially correlated and/or 
non-zero mean. The DMT is shown to depend not on a particular 
fading distribution, but rather on its polynomial behavior near 
zero, and is the same for the simple "amplify-and-forward" 
protocol and more complicated "decode-and-forward" one with 
capacity achieving codes, i.e. the full processing capability at the 
relay does not help to improve the DMT. There is however a 
significant difference between the SNR-asymptotic DMT and the 
flnite-SNR outage performance: while the former is not improved 
by using an extra antenna on either side, the latter can be 
significantly improved and, in particular, an extra antenna can 
be traded-off for a full processing capability at the relay. The 
results are extended to the multi-relay channels with selection 
relaying and typical outage events are identified. 

Index Terms — MIMO, relay channel, outage probability, chan- 
nel capacity, diversity-multiplexing tradeoff. 

I. Introduction 

COOPERATIVE communication strategies have recently 
emerged as an efficient way to exploit multi-user diver- 
sity available in wireless networks to significantly improve 
their performance [l]-[4]. Starting from the pioneering work 
in [1], a number of efficient protocols, distributed space-time 
codes and signal processing strategies have been proposed 
[5]-[16]. While the research was initially concentrated on the 
single-antenna systems [1]-[3][5][12]-[14], the emphasis has 
recently shifted towards multi-antenna systems [6]-[ll],[15]. 
Typical performance metrics in a fading channel include 
the outage capacity or the outage probability [1][11]-[14], 
mean (ergodic) capacity [9] and error rates when specific 
codes/modulation formats are studied [3][10][15]. Due to the 
complexity of the analysis, the performance studies have been 
mainly limited to independent (but not necessarily identically 
distributed) Rayleigh-fading channels. 
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Since the MIMO systems present an additional level of 
difficulty in terms of performance evaluation for all the 
metrics, an elegant framework termed "diversity-multiplexing 
tradeoff (DMT) has been proposed in [17], which allows 
one to quantify the system performance in terms of two 
principal gains, diversity and multiplexing, available in a 
slow-fading MIMO channel when SNR oo [18]. Many 
systems, for which the outage probability/capacity analysis is 
illusive, can be characterized and compared via their respective 
DMTs. This framework has been successfully applied to relay 
channels as well [5]-[8][ll]. Since the amplify-and-forward 
(AF) and decode-and-forward (DP) protocols proposed in [1] 
are not DMT-optimal at high multiplexing gains, a number 
of new protocols were proposed in [5] and their DMT has 
been investigated. This work was further extended to the multi- 
antenna terminals in [6]-[8][ll][34]. While most of the DMT- 
based studies are limited to the SNR oo regime, the finite- 
SNR DMT of the relay channel with independent, Rayleigh- 
fading links and single-antenna terminals have been studied in 
[28]. 

The outage capacity of relay channels with independent 
Rayleigh fading links has been studied in [13] in the low- 
SNR regime, and it was shown that the AF and DF protocols 
are sub-optimal, and a new strategy, the bursty amplify-and- 
forward (BAF), was proposed and shown to be optimal, to 
the 1 ^* order, in the low-SNR low-outage regime, which also 
includes the case of multiple relays. These results have been 
extended to the case of relay channels with generic fading 
distribution (but the links are still required to be independent) 
admitting low-outage Taylor expansion in [14], where it was 
shown that the BAF still achieves the outage capacity to the 
l*** order (but not to the second order) in the low-SNR low- 
outage regime. The case of selection relaying (where the best 
relay only is used) has also been studied and the BAF strategy 
has been shown to be optimal in this case as well. Both of 
these works consider single-antenna terminals. 

In all the studies above, a Rayleigh or Rician-fading channel 
model with independent links has been employed. The only 
exception is [14], where the outage probability/capacity has 
been studied for a generic fading distribution but the analysis 
was limited to the low-SNR regime and the links are still 
required to be independent. 

In this paper, we consider the source and the destination of 
the relay channel equipped with multiple antennas and relay 
nodes equipped with a single antenna (e.g. due to complexity 
constraints). We allow the fading to be non-identical, corre- 
lated Rayleigh or of generic distribution (non-Rayleigh/Rice), 
and consider amplify-and-forward and decode-and-forward 
protocols. The contributions of the paper are as follows: 
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• The outage probability and outage capacity of correlated 
and/or non-identical Rayleigh fading are obtained in a closed 
form (Theorems 1 and 2 in Section III), including insightful 
low-outage approximations (Corollaries 1 and 2). This allows 
one to draw important design guidelines and also to establish 
the limitations of SNR-asymptotic DMT-based designs: two 
systems with identical DMT may have vastly different outage 
performance, e.g. while the DMT of the 1 x 1, 2 x 1 and 
1x2 channels is the same, an additional antenna results in 
significantly lower outage probability (in fact, the ratio of 
1x1 and 1x2 channel outage probabilities grows unbounded 
as SNR — > oo). This is equivalent to a significant SNR 
gain (about 10 dB at the outage probability = lO^'^), not 
captured by the DMT framework. An additional antenna can 
also be traded off for the full processing capability at the relay. 
Furthermore, an extra transmit rather than receive antenna 
is preferable, since, unlike the latter, the former eliminates 
the negative effect of relay noise in the low outage regime. 
Therefore, one concludes that the DMT framework is ill-suited 
for some relay channels. A study of the effect of correlation 
demonstrates that, unlike full-rank MIMO channels, it has here 
a negligible impact on the performance unless the channel 
is nearly-fully correlated. Under certain conditions, the relay 
channel is shown to be equivalent to the maximum ratio 
combining one with an extra array gain. 

• The SNR-asymptotic DMT of the relay channels in the 
AF and DF modes is obtained for a broad class of fading 
distributions including, as special cases, Rayleigh, Ricean, 
Nakagami, and Weibull, which can be non-identical, non-zero 
mean and/or spatially correlated (see Theorems 3-5 in Section 
IV). The AF and DF systems have the same DMT, which 
depends on the minimum diversity order only. This general- 
izes/extends the known DMT or low-SNR results mentioned 
above. 

• These results are extended to the multi-relay channels with 
selection relaying. Typical outage events in relay channels are 
identified. 

The rest of the paper is organized as follows. Section II 
introduces the basic system model. Outage probability and 
outage capacity are studied analytically in Section III. The 
SNR-asymptotic DMT is considered in Section IV for a broad 
class of fading distributions in the AF and DF modes and its 
limitations are pointed out. Multi-relay channels with selection 
relaying are considered in Section V. Finally, Section VI 
concludes the paper and the appendix details derivations. 

II. System Model 

We consider a MIMO relay channel where the source 
(transmitter) and the destination (receiver) are equipped with 
multiple antennas and a relay node equipped with a single 
antenna (see Fig. 1). While both the amplify-and-forward 
and decode-and-forward protocols are considered, the former 
will be assumed for simplicity of exposition, unless indicated 
otherwise, with a fixed gain relay (this is motivated by the fact 
that it is simpler to implement). We assume no direct source- 
destination link. This is motivated by the fact that the direct 
link may be much weaker than the relay one (e.g. no line-of- 
sight etc.) and thus can be neglected (this is the case when 
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Fig. 1 . Relay channel with a single relay node (single-antenna) and multiple- 
antenna source (transmitter) and destination (receiver). 



the relay link is needed most) [3]; it also makes the analysis 
tractable. The case of multiple relay nodes is considered in 
Section V. 

The standard baseband system model of a frequency-flat 
block-fading (quasi-static) relay channel in the amplify-and- 
forward mode is given by. 



y = V KrGrdGsrKdh-tr^ + \/ KrGrdKd£.r + ^ (1) 

where x and y are the source (transmitter) and destination 
(receiver) symbol vectors, h+ and hrd are the source -relay 
and relay-destination normalized channels (i.e. they include 
the multipath fading but not the average path loss), + denotes 
Hermitian conjugation, Gsr and Grd are the source-relay and 
relay-destination average path loss factors, Kr is a fixed relay 
gain, CM{0,<7^) and | - CM{0,alI) are relay and 

destination AWGN noise of variance cTj. and CTq respectively, 
and independent of each other While analyzing the outage 
probability, we assume that h+. and hrd are i.i.d. or correlated 
Rayleigh-fading and are independent of each other, and while 
analyzing the SNR-asymptotic DMT, h+ and hrd are also 
assumed to be independent of each other but with a generic 
fading distribution of a polynomial near-zero behavior (most 
known models are in this class); m and n will denote the 
number of source and destination antennas. Note that the first 
term in ([U represents the signal received at the destination; 
2nd ^jjj yd j-gj-jj^g represent the relay noise propagated to the 
destination and the destination noise. The sufficient statistics 
for y is (see e.g. [18]), 



(2) 



\Kd\ 



= ^KrGrdGsr \Kd\ h+ X + ^ KrGrd \Kd\ + TT^I 

where |h|^ = h+h, and the instantaneous SNR at the desti- 
nation can be expressed as 



(7^ + KrGrd\hrd\'^ CTr ~ l+a\hrd\^ 



7 



7 (3) 



where — xx+ is the covariance matrix of the transmitted 
signal, a = KrGrdCr /cq is the ratio of the average relay noise 
propagated to the destination to the destination noise, and 7 = 
KrGrdGsrO'^/'^o IS the average SNR at the destination, cr^ = 
trUx = x+x is the total transmitted power (at the source). 
The inequality in (O follows from h+ R^jhsr < \hsr\^ cr^, and 
the equality is achieved when R^; — a'^hgrhfr / \hsr\ , i-c. 
beamforming from the source to the relay, x — shgr/ \hsr\, 
where ,s is the scalar transmitted symbol of the total power 
cr^. This, of course, requires channel state information (CSI) 
at the source. When no such information is available, a sensible 
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transmission strategy is isotropic [18], i.e. R^^ — a^l/m. In 
this case, the SNR is 



7 



(4) 



Comparing this to the upper bound in (|3]l, one concludes that 
the source CSl brings in an m-fold SNR gain, but does not 
change the statistics of the instantaneous SNR otherwise and, 
therefore, the outage probability or outage capacity differ by 
a constant SNR shift, and the diversity-multiplexing tradeoff 
is the same in both cases. Below we assume the source CSI 
so that the instantaneous SNR is 



7 



l + a \hrd 
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(5) 



The instantaneous channel capacity (in nats/s/Hz) can now 
be expressed a|3 C — In (1 + 7') and the channel outage 
probability @, i.e. the probability that the channel is not able 
to support a target rate R, is 



Pr{C<i?} = Pr' 



< 



1 



7 



(6) 



l + a |h^d| 

It follows from (|3]l, (|5]l, (|6]l that the optimum (i.e. capacity- 
achieving) transmission strategy is the beamforming at the 
source (towards the relay) and the maximum ratio combining 
at the destination. 

III. Outage Probability, Capacity and DMT at 
Finite SNR 

Since the finite-SNR analysis is not feasible for a generic 
fading distribution, in this section we consider Rayleigh-fading 
links (e.g. source-relay and relay-destination), which may be 
non-identical and/or correlated. First, the outage probability 
is investigated for fixed R. Then, the finite-SNR diversity- 
multiplexing tradeoff is discussed. 

A. Outage probability 

Theorem 1: Let h.rd and h.sr be mutually independent 
circular-symmetric correlated Gaussian random vectors. When 
the eigenvalues of the source-relay and relay-destination link 



correlation matrices R 



sr{rd) 



hsr(rd)h^r(rd 



are non-zero 



and distinct, the outage probability of the single-relay, corre- 
lated channel in the amplify-and-forward mode is 



i=i j=i 



-ax/Xj 



' Ax 



I Ax 



(7) 



where Kn{x) is the iV-th order modified Bessel function of 
the second kind [27], Ai and Bj are the coefficients of a partial 
fraction decomposition given by Ak — Yli^k ^'^/ i^k — Ai) 
and likewise for B, Xi and j]j are the eigenvalues of 'Eisr{rd), 
and X = (e^ — l)/7. 

' Note that the channels in {T) and (2) have the same capacity (since z is 
a sufficient statistics). 

-It is also the best achievable codeword error probability II35II361 . i.e. a 
fundamental performance limit. Realistic codes can be handled via the SNR 
gap to capacity [32]. 



Proof: see Appendix. ■ 

While this expression holds when eigenvalues are different 
and non-zero, the case of identical eigenvalues can be handled 
by the limiting transition using L'Hopital's rule (since the 
outage probability is a continuous function of eigenvalues), 
and zero eigenvalues should be simply excluded, i.e. m, n 
should be treated as the ranks of corresponding correlation 
matrices. 

To get some insight, let us now consider the low-outage 
regime. 

Corollary 1.1: The behavior of Pout in © in the low- 
outage regime x — > (H is determined by the min(TO, n) term 
and is given by: 



Pout ~ aix™ + o(x'"), m < n 
— a2x" + o(a;"), m > n 
= (a3 + 63lni)x'" + o(a;'"), m = n 

where ai, 02, 03 and 63 are constants independent of x, 



(8) 



ai 



02 



as = 



(-1)"+! ^ In 77, 



m! det R. 



det R, 



EE- 

k=l 1=1 



-Dkmia), 



!(n-l)!dctR,d^ A" 

7=1 J 



(-1) 



m+l 



mIdetRsr ' m!(m - 1)! dot R,.d ^ A™ 



In A," 



+63* 



det R, 



b. = 



1 



k=l 1=1 
1 



'-P>km{a), 



m\{m — 1)! det R^^ det Urd ' 



Dkiia) 



(-1) 



l~k^l-k 



[l - k)\{k - l)\kV 

and = ij{k)+ij{k+ll ^dH) = -C, i>{k) - -C+Y!IZI i 
for fc > 2, and C ~ 0.577 is Euler's constant. 

Proof: see Appendix. ■ 

It follows from dHJ that the weakest link (with the lowest 
diversity order) dominates the outage performance, so that the 
channel diversity order is min(m, n). When m = n, both 
links contribute almost equally to the outage and the unusual 
term In i emerges, which has a profound negative impact on 
the outage probability in the low outage regime (a; ^ 1) 
and cannot be found in the regular (no relay) Rayleigh-fading 
channels. Note also that Pout in dS) does not depend on a 
when m > n, i.e. extra Tx antenna(s) eliminate the negative 
effect of relay noise in the low-outage regime. 

In a typical wireless system, the average path loss is about 
50. . . 100 dB or more (unless the transmitter and the receiver 
are veiy close to each other) [23], i.e. Grd ~ 10 
that 



-5 



.10- 



so 



a = KrGrd<yl/(7l < 1 



(9) 



' this requirement means that 7 ^ — 1, which is equivalent to 7 — > 00 
under fixed R (i.e. high average SNR) but also may hold at low SNR as well, 
when < 7 < 1. 
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when the relay gain Kr and <t^Io\ are not too largeQ. 
Motivated by this, we note that the channel in ([T]) and the 
outage probability in (IT) are the same as those of the single- 
keyhole model [24], [25, Theorem 3.1] when a — !• 0, and 
the corresponding keyhole (double-scattering) channel based 
results also apply to the single-relay channel, e.g. the space- 
time codes [26] or diversity-multiplexing tradeoff [19] [20]. 
Under the condition in (|9]l, O can be expressed, after some 
straightforward but lengthy manipulations, as follows 

- A:!(fc - 1)! ^ (A,77,)^ V ^ 

fe=min(m,n) ^ ' ^ ^ ' ' ^ 

(10) 



In the low outage regime, x 
dominates and ( fTOl i reduces to 

^out 



ml{m — 1)! detRi-ddetR, 



m!(m - 



1)! ^ 



0, the min(m, n) term 

-+o(a;'"), TO = n, (11) 
+ o(a;'"), n > m 



detR 



where 6™ = 1/m + 2-(/'([I]), and the m > n case is obtained 
from the n > m one via m ^ n. Clearly, the diversity gain of 
the channel, at fixed rate R, is d = min(m, n). It is also clear 
that the outage probability increases with correlation (since 
det R is maximum for uncorrelated channel and decreases 
with correlation) and the same conclusion holds for (|8). It 
can be shown that the relay noise is negligible in (|8]l and the 
latter reduces to (fTTt when a ^ (m! det Rsr)^^™ and either 
a < {m{m - (for m > 2) or a < J2i=i ^T^m/Vi 

(for m — 1). Note that (to! dot R^^)^^™ decreases with 
correlation, so that a given relay noise may be negUgible when 
the correlation is low, but not when it is high, which is another 
manifestation of the detrimental effect of the correlation. Fig. 2 
shows the outage probability of a 2x2 correlated relay channel 
for different values of a. Clearly, the effect of the relay noise 
is negligible unless a > 1. 

While the results above apply to a correlated channel with 
distinct eigenvalues and do not include the i.i.d. channel, the 
latter can be characterized as follows. 

Theorem 2: Let h^d and hsj.be i.i.d. mutually-independent 
circular-symmetric Gaussian random vectors. The outage prob- 
ability of this single-relay channel in the amplify-and-forward 
mode is 

" ^ " (n-1)! ^ i\{k-i)\ 



-Kn+i-kW ix) 

(12) 

Proof: see Appendix. ■ 

Note that ( fT2b is easy to evaluate numerically since it 
contains finite sums and well-known special functions. It can 
be further expanded as a series in x: 

Corollary 2.1: Pout in Theorem 2 can be expressed as 



Pout = ^ (/;(«) + 9i{oi) Inx) 



^+min(m,n) 



(13) 



i=0 



* As a practical example, Kr = 60 dB, G^d = —100 dB, a"^ = ctq = 
-102 dBm in a typical 3GPP UMTS scenaiio [33], so that a = lO"-* <C 1. 
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Fig. 2. Outage probability of the 2x2 channel vs. a for x = 
10^^, 10~^, 10~^; the normalized Tx/Rx correlation is p = 0.5. Note that 
a < 1 has negligible impact on the outage probability and there is a significant 
increase in the Pout when a > 1, i.e. relay noise is important only when 
it is really high: c^/ctq > Xj [KrGrd]- In a practically-important case of 
= (Tq [33], this reduces to KrG^d > 1 and the relay noise is negligible 
otherwise. 



where fi (a) and gi (a) are independent of x and are given in 
Appendix. The behavior of Pout in the small outage region 
a; — > is determined by the min(TO, n) term: 



Pout 



a^{n — TO + — 1)! 



1)! ^ (TO-fc)!fc! 



+ o(x™), m <n, 



x"{m — n — 1)1 
n!(TO — 1)! 



X" 



+ o(a;"), TO > n, 
a'' In i + b„ 



to! 



(14) 



+ o(x'"), 



(to — 1)! {m — k)lk 

m = n. 

Proof: see Appendix. ■ 

Clearly, the diversity order d — min(TO, n) does not depend 
on a and the unusual term In i is also present, which has 
a profound negative impact on the outage probability. Pout 
does not depend on a when m > n, i.e. extra Tx antenna(s) 
eliminate the negative effect of relay noise in this case as well, 
and increases with a (i.e. with relay noise) otherwise. 

When a — > 0, Pout becomes symmetrical with respect to 
TO and n, which is explained by the symmetry of the channel 
in (HJ in this case, and (fT2] | simplifies to 



1 



in-1) 



1)\ ^ 



k=0 



kl 



1 71 — m I — 1 CO 



4=0 



4=0 



Kn-kiVi^) (15) 

,i+max(m,n) (lna;_c,) 



where /x^. Pi and Ci are independent of x, 

{~iy{\n - to| - i - 1)! 



i!(min(TO,, n) + i)(n — 1)!(to — 1)! ' 

z!(max(TO, n) + i)(\n — m\ + iy.(n — l)!(m — 1)! 

c, = : + + 1) + tp{\n -m\+i + l) 

max(TO, nj + 1 



(16) 
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and ( fT4b reduces tc|f| 



In- 



n > m, 



(17) 



m!(m — 1)! 

=-f-y 

m! VG/ 

where G — {{n — l)...{n — m))^/'" > 1 and the m > n case 
is obtained from the n > rn one via the substitution m ^ n. 
Note that Po„t = Pmrc{x/G) < Pmrc{x) when n > m, 
where Pmr.c{x) ~ x'"/m! is the outage probability of m- 
branch maximum ratio combiner (MRC) in the i.i.d. Rayleigh- 
fading channel, so that the relay channel is better than (n > 3) 
or equal to (n — 2) the m-branch MRC channel in this case 
and, as far as the outage probability is concerned, the relay 
channel is equivalent to a cascade of a fading link (source- 
relay) and a non-fading relay-destination Unk with an SNR 
gain of beamforming at the destination equal to G. When n ^ 
m, G ~ n, i.e. a gain of n-element antenna array. Similar 
conclusions hold for the n < m case via the substitution m -H> 
n. 

B. Special Cases 

To obtain some insight, let us now consider 1 x 1, 2 x 1 
(2 Tx, 1 Rx antenna) and 1x2 i.i.d. channels. In these cases, 
(fT2] i simplifies to 

nlxl 



P, 



out 

)2xl 

out 

1x2 



e 

2xe-"^ 



4xi^i(v4x), 

K2{V4^). (18) 



At the low outage regime ax <C 1, this can be approximated 
as 



Pi 



x{a + \n^), Pif «a:(l + a), P^V 



X, 



(19) 



so that Pi,r/ > Pif 



> p2xl nnH plxl/pl> 



Pout^/Pout^ grow unbounded as a: 0, even though the 
diversity gain = 1 in all three cases. Note also that P^^/ = 
Pr < Pout^ in the low outage regime, where Pr w x 
is the outage probability of the Rayleigh 1x1 channel 
(no relay), i.e. the 2x1 relay channel is equivalent to the 
classical 1x1 Rayleigh channel (no relay) but the 1x2 
one is not, unless the impact of relay noise is negligible (i.e. 
a = KrGrdCr /o'o ^ !)■ This is the case in ( fT9] l when 
a ^ In i and a ^ 1 for the 1x1 and 1x2 channels 
respectively, and it is always negligible for the 2x1 channel in 
the low-outage regime (we attribute this to the higher diversity 
order of the source-relay link in this channel). 

As Fig. 3 demonstrates, the approximations in (fT9] l are 
accurate for a; < 1 (low outage regime). They also provide 
an insight into the way typical outage events occur in the 
relay channels: for 1x2 and 2x1 channels, it is when the 
lowest diversity order link is in outage (S-R and R-D links 
respectively) and the other link's contribution is negligible. For 

^ It can be shown that the impact of relay noise is negligible so that (TT) 
can be used if a ^ 1/ (m(n — m)) for n > m (the normalized contribution 



of the relay noise does not exceed m{n 
for n = m. 



■ m)a in this case), and if a <C 1 




Fig. 3. Outage probability of 1 X 1, 2 X 1 and 1x2 relay channels in 
(18) with corresponding approximations in l |19t for a = 1. The low-outage 
approximations are accurate for x < 1, and Monte-Carlo simulations validate 
the analysis. 



1x1 channel, it is not the case anymore since P^^/ were about 
(2 + a)x otherwise. Furthermore, (2 + a)x <C x{a -f- In i) as 
X — > 0, which tells us that the typical outage events are when 
both links are in partial outage. 

Let us now compare the m x m relay channel to the to x 1 
MRC channel when a — > 0, which has the same diversity 
order and whose outage probability Pmrc ~ x"^ /ml at the 
low outage regime, so that Pout Pmrc ~ ^'n^/im - 1)!, 
which also grows unbounded at low outage /large SNR regime, 
even though they have the same diversity order Thus, the tox 1 
MRC channel performs much better than to x to single-relay 
channel at low outage. Comparing the MRC channel to (to + 
1) X TO relay channel, Pi,™^^''^"V^A/JJC ~ 1/to! < 1, i.e. the 
(to + 1) X TO relay channel performs better, which re-enforces 
our earlier conclusion that one extra source/destination antenna 
improves the performance significantly, even though it does 
not affect the diversity order (and also the DMT - see Section 
IV). Typical outage events can be identified in a similar way. 

To explore the impact of correlation, let us now consider 
the 2x1 correlated relay channel when a — > 0. In this case 



1 p 



where p is the (scalar) normalized correlation. The eigenvalues 
'71,2 = 1 ± and the outage probability in ^ can be 
explicitly expressed as 



Po 




(20) 



where the approximation holds in the low outage regime. 
While the effect of correlation is detrimental in general, it 
is significant only when the correlation is very high - see Fig. 
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Fig. 4. Outage probability vs. normalized correlation for 2x1 AF single- 
relay, Rayleigh-fading channel (a = 0), at z = 10~^ (at R = Ibit/s/Hz, this 
con'csponds to 7 = 20di3). Note that coiTelation has significant effect only 
when it is very high, \p\ > 0.8... 0.9, and that the approximation in )20t is 
accurate unless |p| — > 1. Compared to 1x1 channel, the additional destination 
antenna significantly reduces the outage probability, even for high correlation. 
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Fig. 5. Outage probability vs. SNR for 1x1, 2x1 and 1x2 AF single-relay 
i.i.d. Rayleigh-fading channels; r = 0, a = 1, 10® Monte-Carlo trials. While 
the SNR-asymptotic DMT is the same for all channels, the gap between the 
outage probabiHties of the 1x1 and 2x1 (1x2) channels increases with SNR; 
at SNR=40dB, the difference is about lOdB. Unlike an extra Rx antenna, an 
extra Tx one allows to combat relay noise. 



4 , and the outage probability increases only logarithmically 
with correlation. Comparing to the 1x1 channel, an additional 
Rx antenna significantly reduces the outage probability, even 
when the correlation is high. When p ^ 0, ( |20] | reduces to 
Pint in (O. as it should. 

C. Finite-SNR DMT 

At finite SNR, the definitions of the diversity and multiplex- 
ing gains in [17] are used without the limiting transition (see 
e.g. [28][29]), 

d= -lnFo„t/ln7, r = i?/ 1117 (21) 

from which it follows that x — {e^^ — ~ 1/7^^'' and 
the approximation holds at sufficiently high (but finite) SNR, 
7^"'' :$> 1. Using this in (fT4] i. one obtains for 1 x 1, 2 x 1 and 
1x2 i.i.d. channels and < r < 1 : 

PiV«(« + (l-01n7)/7'-^ 

Pif « (1 + a)/7l-^ PiV « Ih'-"- (22) 

Thus, the diversity gain can be expressed at high SNR/low 
outage as 

c^ixi ~ 1 ^ J' — In (a + (1 — r) ln7) /ln7 
dix2 w 1 - r — ln(l + a)/ ln7 
^2x1 w 1 - r 

It is straightforward to see that the SNR-asymptotic DMT (i.e. 
when 7 — > 00) is the same in all three cases: d{r) — 1 — r. 
However, the outage probabilities behave quite differently, as 
Fig. 5 demonstrates. From (l22t . 

|f|«a+(l-r)ln7, ^ ^ I + a, (23) 

out out 

so that PiV/^ol?, ^ontV^o«f ^ 00 as 7 ^ 00 for 
r < 1, i.e. there is a significant advantage in using an extra 
antenna at either Tx or Rx end at high SNR, even though 



the SNR-asymptotic DMT is the same in all three cases. 
Thus, the corresponding DMT-based conclusion breaks down 
in a most dramatic way. The ultimate reason for this is the 
In 7 term in ( |22] |. not captured by the SNR-asymptotic DMT, 
which indicates that the latter is ill-suited for relay channels. 
Furthermore, (|23) also shows that an extra Tx rather than 
Rx antenna is preferable since, unlike the latter, the former 
eliminates the effect of relay noise at the low outage regime. 

When the relay node has full processing capability, i.e. the 
decode-and-forward protocol with capacity-achieving codes on 
both ends, the relay channel capacity is (7 = min {Csr, C'rd}, 
where Csr,Crd are the capacities of the source -relay and 
relay-destination links, so that the weakest link dominates the 
outage performance, and 

plxl ^ 2 1x2 _ p2xl ^ ^ 

^ out ~ i^r ' out ^ out ~ i-r \^^> 

assuming for simplicity the same average SNR on the source- 
relay and relay-destination links. Clearly, the SNR-asymptotic 
DMT d{r) — 1 — r is the same in all cases, and also the same 
as for the AF mode (this is further generalized in Section IV to 
a broad class of fading distributions). Comparing (l24t to (l22b . 
we note that the 2x1 channel in the AF mode outperforms 
the 1x1 channel in the DF mode, i.e. an extra Tx antenna 
in the AF mode is better then the full processing capability at 
the relay and can be used as a simple alternative of the latter, 
even though it does not improve the SNR-asymptotic DMT. 
The same can be said about an extra Rx antenna if a < 1. 
This re-affirms the earlier conclusion that the SNR-asymptotic 
DMT framework is ill-suited for relay channels and should be 
used with extreme caution when formulating design guidelines 
and designing space-time codes. It also follows from (|24] | that 
an extra antenna in the DF mode brings a modest SNR gain 
(= 3dB at r = 0), unlike the AF mode where this gain can be 
very significant (see Fig. 5). Comparing (l24t to (l22t . we also 
note that the full processing capability at the relay does not 
bring in any advantage for the 2x1 channel; the same applies 
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to the 1x2 channel when a ^ 1. Clearly, a comparison of 
different channels/systems in terms of the outage probability 
may agree with the SNR-asymptotic DMT-based one in some 
cases while significantly disagreeing in others. 



D. Outage Capacity 

Based on the outage probability results of the previous 
section, we can now analyse the outage capacity which is 
defined as = max {R : Pout{R) < e} and can be expressed 
as [18] 

Ce =ln(l + 72;^) 

Ce ~ ln7 - In — , > 1 (high SNR), (25) 

Ce W JXe, JXe < 1 (low SNR), 

where = Pg^itC^) is the inverse of Pout{x), which quantifies 
the SNR loss compared to the AWGN channel whose capacity 
is Cawgn = In (1 + 7). Following ([8]l, the SNR loss can be 
approximated as 



(e/ai) 
(e/aa) 





m < n 




m > n 






e ■ m 





(26) 



03 



6, In ^ 



i.e. it scales roughly as e^/ mm(m,n) ^jjj^ jj^g outage probability 
and increasing min(m, n) has a significant positive effect on 
the outage capacity, especially in the small outage region. 
Since the non-fading AWGN capacity is Cawgn ~ ln7 (high 
SNR) and Cawgn ~ 7 (low SNR), the capacity loss in the 
fading relay link for m > n is 



AC — Ce — Cawgn 



-In— (high SNR) 

n e 

A = C, /Cawgn = (£/a2)'/" (low SNR) (27) 



i.e. an additive loss at high SNR and multiplicative at low, so 
that this effect is much more severe in the latter case. Similar 
conclusions can also be obtained for the n > m and n ~ m 
cases. 

To obtain some insight, let us compare 1x1 and uncorre- 
cted 1 X 2, 2 X 1 channels 



1 ' 



1 X 1 



a + In ■ 

■ e/{l + a), 1x2 
2x1 



(28) 



The ratio of the SNR loss factors for 2 x 1(1 x 2) and 1x1 
channels grows unbounded as In i when e — > 0, i.e an extra 
Rx or Tx antenna has a significant positive impact. Fig. 6 
compares the normalized outage capacities of 1 x 1 and 2x1 
channels, clearly indicating a significant benefit of an extra 
antenna, especially in the low SNR regime. 




30 40 
SNR [dB] 

Fig. 6. The outage capacity of 1x1 and 2x1 channels normalized to the 
AWGN channel capacity versus SNR; solid line - Monte-Carlo simulations, 
dashed line - via the approximation in 1261 and l l25t . While an extra Tx 
antenna brings 70% increase in the outage capacity at high SNR = 40 dB, it 
brings much more dramatic, 5-fold increase at low SNR = -10 dB. 



IV. SNR-AsYMPTOTic DMT 

Following [17] [18] [21], we define the SNR-asymptotic di- 
versity and multiplexing gains as 



lim lnPo„t/ln7, r= Km i?/ln7 (29) 

7— >oo 7— foo 



Using the results of Section III, the SNR-asymptotic diversity- 
multiplexing tradeoff in the single-relay correlated channel can 
now be characterised. 

Theorem 3: Consider the relay channel in ^ under the 
conditions of Theorem 1 (correlated non-identically distributed 
fading) or Theorem 2 (i.i.d. fading). Its SNR-asymptotic DMT 
can be expressed as 



rf(r) = min(TO, n)(l — r), 0<r<l 



(30) 



Proof: it follows from (l29l l that R ~ rlnj and x = 
(e^ — l)/7 = 1/7^"*^ as 7 — >■ 00. Using this in (|8]l to obtain 
Pout and substituting it in ( l29l) . one obtains d as in ( l30l l. In 
the case of i.i.d. fading, ( fT4l ) is used instead of ([8]) and the 
same result is obtained. ■ 

Theorem 1 in [21] (for single relay case) is a special case 
of Theorem 3 here, i.e. when hgr and h^d are i.i.d. complex 
Gaussian. The DMT of relay channels/networks under i.i.d. 
Rayleigh/Rician fading has been also studied in [34] [35]. 
Following a different line of analysis [30] [31], a significant 
generalization of this result can be obtained for a broad class 
of fading distributions (including, as special cases, Rayleigh, 
Rice, Nakagami, and Weibull, which may be correlated, non- 
zero mean and non-identically distribute^!). 

Theorem 4: Consider the relay channel in ([TJ such that the 
PDFs js(x) and of |hsr|'^ and |hrci|^ behave polynomi- 

ally near zero, i.e. ]s{x) ~ x'^^'^^ , fd{x) ~ x'^''~^ as x — > 0, 
where ds, dd are the diversity gains (orders) of the source-relay 
and relay-destination links at r = 0, and f{x) ^ g{x) means 

' It is straightforward to show that full-rank correlation does not affect the 
order of the polynomial behavior of | h^r | or | h^^ | near zero. 
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that there exist positive constants < A < B < oo, such that 
Agi^) < f{x) < Bg{x) for sufficiently small x. The DMT 
of this channel in the amplify-and-forward mode is 

d(r)=min(4(r),dd(r)) (31) 
= min(4,dd)(l -r), 0<r<l 

where ds[r) = ds{l — r),dd{r) — d^il — r) ai-e the SNR- 
asymptotic DMTs of the source-relay and relay-destination 
links. ■ 

Theorem 4 demonstrates that the SNR-asymptotic DMT 
is affected by the number of degrees of freedom (diversity 
order) available in the channel and not by particular fading 
distribution, as long as the definition of diversity gain in (|29] l 
makes sense. Similar result has been also established for full- 
rank MIMO channels in [22]. It follows from Theorem 4 that 
full-rank Tx/Rx correlation and the relay noise do not affect 
the SNR-asymptotic DMlQ Furthermore, from (g) and (|5), 
the absence of CSI at the source is equivalent to an m-fold 
SNR loss and, therefore, has no effect on the SNR-asymptotic 
DMT (i.e. ( |3TI ) still holds). The transmit beamforming in 
combination with QAM modulation (e.g. see [18]) is an 
example of a space-time code that achieves the DMT of the 
single-relay channel with the source CSI. When no such CSI 
is available, isotropic transmission in combination with QAM 
will achieve the DMT. 

Let us now consider the DMT of the decode-and-forward 
single-relay channel, assuming capacity achieving codes and 
complete decoding/encoding at the relay. Following similar 
line of analysis, one can establish its DMT for the broad class 
of fading distributions. 

Theorem 5: Under the conditions of Theorem 4, the 
diversity-multiplexing tradeoff of the single-relay channel in 
the DF mode is the same as in the AF one and is given by ( [3T1 ), 
i.e. full processing capability at the relay does not improve the 
DMT. ■ 

Thus, the single-relay channel subject to fading from a broad 
class of distributions has the same SNR-asymptotic DMT 
in the amplify-and-forward and decode-and-forward modes. 
However, as we argued in section III, the full processing 
capability at the relay does help to reduce significantly the 
outage probability when the SNR is finite. Furthermore, fol- 
lowing Theorem 3, the DMT is affected by min(TO, ?i so that 
from the SNR-asymptotic DMT viewpoint there is no point in 
using unequal number of antennas (e.g. no point to use more 
source than destination antennas). This, however, does not 
hold when SNR is finite as our analysis in Section III shows: 
an additional source or destination antenna can be traded-off 
for the full processing capability at the relay and improves 
the performance significantly. Therefore, the SNR-asymptotic 
DMT should be used with caution when formulating design 
guidelines. This conclusion also applies to full-rank MIMO 
channels for a broad class of fading distributions [29]. 

^ This has been established for Rayleigh-fading keyhole channels in [24]. 
* In the case of i.i.d. Rayleigh and Rician-fading links, this conclusion was 
obtained in [21] and [35]. 



V. Selection Relaying 

Let us now consider the case of multiple relay nodes 
and when the selection relaying is used, i.e. only the best 
relay node is used at any time, out of N available nodes. 
This protocol was originally proposed in [16] and was also 
considered in [14] in the low-SNR regime. It is motivated by 
its low complexity and also by the fact that little interference is 
created to other users since only one relay is transmitting. We 
assume that different relay links are independent of each other 
(which is justified by geographical separation of the relays), 
that each link follows the statistics of section III (for outage 
probability analysis), and that amplify-and-forward protocol is 
usecH (unless indicated otherwise). 

A. Outage Probability 

An outage takes place when the best relay link and therefore 
all the relay links are in outage, so that the outage probability 
is 

N 

Pout=Pr{C<i?} = []P, (32) 

where Pi = PrjC^ < R] is the outage probability of z-th 
relay link, which is given in (I7]i-(l20ll. and each link is allowed 
to have its own statistics. When all the links have the same 
statistics, ( l32l ) reduces to 

Pout = Pi (33) 

so that the outage probability decreases exponentially in N 
compared to the single relay case, which is especially signif- 
icant in the low outage regime, when Pi <C 1. From (l33l l 
and (HI), (O, the diversity order \s d — N min(m, n) under 
Rayleigh fading (also with full-rank correlation), so that the 
number of antennas and relay nodes can be traded off for each 
other. 

To obtain further insight, let us consider 1 x 1, 1 x 2 and 
2x1 i.i.d. channels, for which ( |33] | simplifies in the low outage 
regime to 

Plx/«x^(a + lni)^, (34) 
Pout ~ ^ (1 + a) , Pout ~ ^ 1 

so that the ratio P^ut^ / Pint"" ~ + grows un- 

bounded as a; — T' (i.e. 7 -> oo under fixed P or P — > 
under fixed 7) and, as in the case of iV — 1, there is a 
significant advantage in using 2 source antennas instead of 
1 at this regime, even though the diversity gain is the same in 
both cases. This conclusion also holds when the source and 
destination are equipped with more antennas. 

B. Diversity-Multiplexing Tradeojf 

Based on the outage probability in (l32T i and Theorems 
3, 4, the SNR-asymptotic DMT of selection relaying can 
be immediately characterized for a broad class of fading 
distributions. 

' In this case, the best relay is arg maxi | Ih^jj j |^ Ihg^.i P /(I + |h„j_i |^) 
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Theorem 6: Consider the relay channel in ([TJ under the 
conditions of Theorem 4 and 5 and assume that all relay 
links are independent of each other Its diversity-multiplexing 
tradeoff in the AF and DF modes is as follows: 

N 

d{r) = ^ mm{ds,i{r), dd,i{r)) (35) 
1=1 

N 

= (1 - r) ^ min(ds_i, < r < 1, 

i=l 

where ds,i{r) = - r),dd,t{r) = dd,i{l - r) are the 

DMTs of the Tx-relay and relay-Rx links. ■ 

It follows that the total DMT is the sum of the DMTs 
for each relay and increasing the maximum of {ds^i,dd.i) (by 
increasing the number of antennas at the corresponding end 
of the link) will not improve the DMT if mm{ds,i, dd.i) are 
fixed, which is the same as for the single relay case. 

In the case of identical link statistics, ( |35] ) simplifies to 

d{r) = N{l-r)mm{ds,dd), 0<r<l (36) 

and A^-fold increase in the diversity gain is obvious. However, 
we caution that the same limitations of DMT-based analy- 
sis/design as for the single-relay case also hold for selection 
relaying. 

Finally, we comment that the impact of direct link (source- 
destination) can also be included in the analysis in the same 
way: any reasonable relaying protocol makes use of both links 
so that an outage takes place when both links are in outage 
and the overall outage probability is PoutPdu where Pdi is the 
outage probability of the direct link (this is exactly the case 
when selection relaying is used), and the overall diversity gain 
is the sum of per-link diversity gains. 

VI. Conclusion 

Outage probability and diversity-multiplexing tradeoff have 
been investigated for MIMO relay channels. The SNR- 
asymptotic DMT has been established for such channels under 
a broad class of fading distributions, thus generalizing earlier 
results known for i.i.d. and correlated Rayleigh channels. The 
latter two have been investigated in greater depth. Compact, 
closed-form expressions, and corresponding low-outage ap- 
proximations have been obtained for their outage probability 
and capacity, which, unlike the SNR-asymptotic DMT, also 
hold at realistic SNR values. Significant difference between 
the SNR-asymptotic DMT and finite-SNR outage performance 
has been emphasized. In particular, while the SNR-asymptotic 
DMT is not improved by using more antennas on either side, 
the outage probability can be significantly improved and, in 
particular, an extra Tx antenna can be traded-off for full 
processing capability at the relay. The results are extended 
to channels with multiple relays under selection relaying. The 
SNR loss of fading relay channels compared to the AWGN 
channel and the impact of relay noise have been quantified. 
Under certain conditions, the relay channel has been shown to 
be equivalent to the maximum ratio combining channel. 
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Appendix 



Let Qs = |hsr|^, Qd = jhrdl^- Since they are independent 



A. Proof of Theorem 1 

Let gs = \h.sr\ 
and non-negative. 



Pmit — Pr 



gs9d 



< X 



1 + agd 

OO px{l-\-OLt-i) jti 

fd{ti) / fs{t2)dt2dti, 



(37) 



where fs{t) and fd{t) are PDFs of gs and g^. Under the 
adopted assumptions, the distribution of g is generaHzed 
with the characteristic function ^i3{uj) = dct [I — jwR]^^, 
where R is the correlation matrix and j = \/^^. When 
R is non-singular and has N distinct eigenvalues Xk, the 
characteristic function (CF) of g is 



N 



N 



B. Proof of Corollary 1.1 

Using the series expansion of Km{x) [27], 

N-l 

2 



1 JV-i 
Kn{x) = -Y.[-1) 



AN-k-l)\ 



k=0 



kl{x/2) 



N-2k 



(40) 



f 1 (2^/2) 



N+2k 



k=0 
X 1 



k\{N + ky. 



1 



X I ln---V;(fc + l)--VXA^ + fc + l) 



and making use of the following properties of partial fraction 
decomposition. 



N 



N 



Y^Ak^l, ^^=0, ^ = l,...,7V-l, (41) 

k=l k=l 

one obtains after lengthy but straightforward manipulations an 
alternative expression for (|7): 
(i) when m > n. 



00 m 



-Tout - 2^ 2^ ^, j^k-^ 

k—7n j — 1 3 

k=nl=k t,] '^j^i ^ 



(42) 



00 00 



E E E 



Aj Bi In rji 



Dki{a)x^ 



k=l i=max(fc,m) i-j 3 * 

(ii) when m < n 



Pout 



k—'rn j — 1 J 



(43) 



00 00 



A,B,, 



■E E Ei3'^-(")-'0-f 

k—l l—ma-ic{k,m) i,j ■? ^ 



00 00 



^ — T V ^ — ^ Aj Bi In. \j 
Z-^ Z-^ Z-^ }J-fj^ 

k—n l—k i,j 3 



Dki{a)x' 



The leading term in the above series corresponds to fc = m 
when m < n, and k = n when m> n. Keeping only this term 
and making use of the following property of partial fraction 
decomposition: 



k=l ^fe 



where Ak are the coefficients of the partial fraction decompo- 
sition of ^g{uj). Thus, the PDF of g is: 



one obtains ^ after some manipulations. 



fc=i 



(39) 



where a; > 0. After substituting ( [39] l into (|37| i and integrating 
using the standard integrals [27], ([7} follows. ■ 



C. Proof of Theorem 2 

Under the adopted assumptions, g is the central random 
variable with the following PDF: fg{x) = x'^-^e''' / {N -1)1, 
where N (to or n) is the number of degrees of freedom. After 
substituting this into dJTl l. integrating (using standard integrals 
in [27]), and making some manipulations, ( fT2] l follows. ■ 
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ci{a 

//'(« 
c;(a 



1=0 j = [i+i_p.]^ 

E E ^«.(") 

j=l j = + 
inin(;— po.m) l—pQ—i 

E E "'^(« 

inin(/— po,TTi) l—pQ—i 



(-l)'(po-gz,, -1)! 



Q (a), m < n, 
ci{a), m = n; 

I I ^ Po, and otherwise; 



(45) 



E E i-^r^'^'^^M) 



[qiij -poY-qiojKi + 

i>{mij - po + 1) + 'iPiqioj + 1) + + 



i\jl{n — l)!(r7i — i)! 



min(i,po-l) l-i 

E E ^^^■(") 

i=0 j = [l+l-po] + 



{quj - poY-qiojKl + m) 

, p. = \m-n-i\, qiij=l-i-j; 
{-iy-^{po-qioj-iy. 



, I ^ Po, and otherwise; 



qiijl{l + n) 



(46) 



Qiji - 



m l~po (—lY^P°(— —1)1 

V V r!.,(a)i-^ ^ ,^'7 , l>Po, and otherwise; 

^=p„+lJ = [l+l-^]+ ^^"J-J ^U'' ^ ^ 

inin(/,m) l — i min(/,m) l—i 

E E'^*J'^*j(")' = E E®y'^*j(")' ^ > Po + 1, and otherwise; 

j=Po + l J=0 i=po + l i=0 

Po '-Po Po '-Po 

E E '^iji^ijic')^ 9i{c() = E E Qiji^ijio^)^ I > Po, and otherwise; 

1=0 j"=0 2=0 j=0 

(__l)p.+i+i 



qiijKqioj -po)!(^ + "-)' 



T,j; = + 1) + V'('7iOi - Po + 1) + + ") ^) ©iji 



D. Proof of Corollary 2.1 

Using the series expansion of Kn{x) (see ( |40] |) and ( fT2b . 
one obtains after lengthy but otherwise straightforward ma- 
nipulations, ( fTsT l, where /((a) and (7i(Q!) are independent of 
X. When m < n fi, (a) and gi (a) are given by (l45T l. where 
[a;]+ = a; if X > and otherwise. When m > n, 

= //(a) + //'(a) - c'lia) - q(a), 



cj(a) and c"(q!) 



are 



where //(a 



//'(a), . 
independent of x, and given by (146b . The first non-zero term 
in the above series corresponds to Z = 0, so that the lowest 
pow er of X is min(m,n). By keeping this term, one obtains 
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